Abstract. We compute the temperature profiles of accretion discs around rapidly ro- 
tating strange stars, using constant gravitational mass equilibrium sequences of these 
objects, considering the full effect of general relativity. Beyond a certain critical value 
of stellar angular momentum (J), we observe the radius (r or b) of the innermost stable 
circular orbit (ISCO) to increase with J (a property seen neither in rotating black holes 
nor in rotating neutron stars). The reason for this is traced to the crucial dependence of 
' dr OT i, I dJ on the rate of change of the radial gradient of the Keplerian angular velocity 
O ; at , orb with respect to J. The structure parameters and temperature profiles obtained 
C^l ■ are compared with those of neutron stars, as an attempt to provide signatures for distin- 
q^' guishing between the two. We show that when the full gamut of strange star equation of 
state models, with varying degrees of stiffness are considered, there exists a substantial 
^ ! overlap in properties of both neutron stars and strange stars. However, applying accretion 
disc model constraints to rule out stiff strange star equation of state models, we notice 

i— I ■ 

■ that neutron stars and strange stars exclusively occupy certain parameter spaces. This 
IT) ' 

qq ' result implies the possibility of distinguishing these objects from each other by sensitive 
I observations through future X-ray detectors. 
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1. Introduction 

Low mass X-ray binaries (LMXBs) are believed to contain either neutron stars (NSs) 
or black holes accreting from an evolved or main sequence dwarf companion that fills 
its Roche-lobe. The proximity of the companion in these systems cause matter to spiral 
in, forming an accretion disc around the central accretor. Observations of LMXBs can 
provide vital clues of the structure parameters of the accretors and, in particular for 
NSs, this can lead to constraining the property of the high density matter composing 
their interiors. Therefore, the estimation of the radius of the central accretor in SAX 
J1808.4 and 4U 1728-34 (Li et al. 1999a; Li et al. 1999b; Burderi & King 1998; Psaltis & 
Chakrabarty 1999) indicating the object to be more compact than stars composed of high 
density nuclear matter, acquires significance. These results moot alternate suggestion 
about the nature of the central accretors in at least some of the LMXBs. 

In this regard, the strange matter hypothesis, formulated by Bodmer (1971) and Wit- 
ten (1984) (see also Itoh 1970; Terazawa 1979), has received much attention recently. 
The hypothesis suggests strange quark matter (SQM, made up of u, d and s quarks), in 
equilibrium with weak interactions, to be the actual ground state of strongly interacting 
Send offprint requests to: Sudip Bhattacharyya 
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matter rather than 56 Fe. If this were true, under appropriate conditions, a phase tran- 
sition within a NS (e.g. Olinto 1987; Cheng & Dai 1996; Bombaci & Datta 2000) could 
convert the entire system instantaneously into a conglomeration of strange matter or, as 
is commonly referred to in literature, strange stars (SSs). 

It is of fundamental interest - both for particle physics and astrophysics - to know 
whether strange quark matter exists. Answering this question requires the ability to 
distinguish between SSs and NSs, both observationally as well as theoretically and this 
has been the motivation of several recent calculations (Xu et al. 2001; Gondek-Rosinska 
et al. 2000; Bombaci et al 2000; Zdunik 2000; Zdunik et al 2000a; Zdunik et al 2000b; 
Datta et al 2000; Stergioulas et al 1999; Gourgoulhon et al 1999; Xu et al 1999; Gondek & 
Zdunik 1999; Bulik et al. 1999; Lu 1998; Madsen 1998). One of the most basic difference 
between SSs and NSs is the mass-radius relationship (Alcock, Farhi & Olinto 1986): while 
for NSs, this is an inverse relationship (radius decreasing for increasing mass), for SSs 
there exist a positive relationship (radius increases with increasing mass) q In addition 
to this difference, due to SSs being self-bound objects, there also exists the possibility 
of having configurations with arbitrarily small masses; NSs on the other hand, have 
a minimum allowed mass (e.g. Shapiro & Teukolsky 1983; Glendenning 1995; and more 
recently, Gondek et al. 1997; Gondek et al. 1998; Goussard et al. 1998; Strobel et al 1999; 
Strobel & Weigel 2001). Nevertheless, it must be remarked that for a value of gravitational 
mass equal to 1.4 M Q (the canonical mass for compact star candidates), the difference 
between the predicted radii of nonrotating configurations of SS and NS amounts, at most, 
only to about 5 km; a value that cannot be directly observed. There arises, therefore, a 
necessity to heavily rely on models of astrophysical phenomena associated with systems 
containing a compact star to estimate the radius: for isolated pulsars, models of glitches 
(e.g. Datta & Alpar 1993; Link et al. 1992) have been used in the past for making 
estimates of the structure parameters and for compact stars in binaries, such estimates 
have been made by appropriately modelling photospheric expansion in X-ray bursts 
(van Paradijs 1979; Goldman 1979) and more recently by constraining the inner-edge of 
accretion discs and demanding that the radius of the compact star be located inside this 
inner-edge (Li et al. 1999a; Li et al. 1999b; Burderi & King 1998; Psaltis & Chakrabarty 
1999). In particular, the work by Li et al. (1999a; 1999b) suggest strange stars as possible 
accretors. However, these calculations did not include the full effect of general relativity. 
Even on inclusion of these effects (Bombaci et al 2000), the results for at least one 
source: 4U 1728-34, remain unchanged. There have also been contradictory reports on 

1 In the present work, we consider only bare strange stars, %. e. we neglect the possible presence 
of a crust of normal (confined) matter above the deconfined quark matter core (see e.g. Alcock, 
Farhi & Olinto 1986). 
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the existence of strange stars: for example, calculations of magnetic field evolution of 
SSs over dynamical timescales, make it difficult to explain the observed magnetic field 
strengths of isolated pulsars (Konar 2000). On the other hand, Xu & Busse (2001) show 
that SSs may possess magnetic fields, having the observed strengths. These magnetic 
fields, these authors argue, originate due to dynamo effects. In our analysis here, we 
ignore the effects of magnetic field. 

In this paper, we calculate constant gravitational mass equilibrium sequences of rotat- 
ing SSs, considering the full effect of general relativity. We solve Einstein field equations 
and the equation for hydrostatic equilibrium simultaneously for different SS equations of 
state (EOS) models. We compare our theoretical results with those obtained previously 
for NSs (Bhattacharyya et al. 2000). In addition, we calculate the radial profiles of ef- 
fective temperature in accretion discs around SSs. These profiles are important inputs 
in accretion disc spectrum calculations, crucially depending on the radius of the inner 
edge of the accretion disc. This radius is determined by the location of r or b with respect 
to that of the surface (R) of the star, both of which are sensitive to the EOS, through 
the rotation of the central object. In particular, we notice that r OID increases with stellar 
angular momentum (J) beyond a certain critical value (a property not seen in either 
rotating black holes or neutron stars). We trace this behaviour to the dependence of 
dr orb I dJ on the rate of change of the radial gradient of the Keplerian angular velocity at 
r or b with respect to J. The prospect of using the temperature profiles for calculation of 
accretion disc spectrum and subsequent comparison with observational data, therefore, 
gives rise to the possibility of constraining SS EOS, and eventually to distinguish between 
SSs and NSs. 

The structure of the paper is as follows. In section 2, we describe briefly our method 
of calculation of the structure of rapidly rotating relativistic stars and the accretion disc 
temperature profile. We also provide here, a brief description of the EOS models used. 
The results described in section 3 are discussed and summarised in section 4. 

2. Methodology 

2.1. Structure calculation 

We calculate the structure of rapidly rotating relativistic stars using an iterative proce- 
dure described in Cook et al. (1994) (see also Datta et al. 2000). The metric describing 
the space-time around a rotating relativistic star can be given as 
dS 2 = g^dx^dx^i^u = 0,1,2,3) 
= -e'+Pdt 2 + e 2a (df 2 +f 2 de 2 ) 

+e^- p f 2 sin 2 6{d<p - ujdt) 2 (1) 
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This metric is stationary, axisymmetric, asymptotically flat and reflection-symmetric 
(about the equatorial plane). The metric potentials -f,p, a and the angular speed (uj) 
of zero-angular-momentum-observer (ZAMO) with respect to infinity, are all functions 
of the quasi-isotropic radial coordinate (f) and polar angle (9). The quantity f and the 
Schwarzschild-like radial coordinate (r) are related through the coordinate transforma- 
tion: r — fe( 7 ~ p )/ 2 . We use the geometric units c = G = 1 in all the equations that 
appear this paper. 

With the assumption that the star is rigidly rotating and described by a perfect 
fluid, we solve Einstein field equations and the equation of hydrostatic equilibrium self- 
consistently and numerically from the centre of the star upto infinity to obtain 7, p, a, uj 
and f2* (angular speed of neutron star with respect to an observer at infinity) as functions 
of f and 9. The inputs for this calculation are: the chosen EOS, assumed values of central 
density and ratio of polar to equatorial radii. The outputs are bulk structure parameters 
such as: gravitational mass (M), equatorial radius (R), angular momentum (J), moment 
of inertia (/) etc. of the compact star. We also calculate the specific disc luminosity (-Ed), 
the specific boundary layer luminosity (Ebl), the radius (r or b) of the innermost stable 
circular orbit (ISCO), specific energy (E), specific angular momentum (l) and angular 
speed (fix) of a test particle in Keplcrian orbits (Thampan & Datta 1998). 

In the present work, for each adopted EOS, we construct constant M equilibrium 
sequences with f2* varying from the non-rotating case (static limit; f2* = 0) upto the 
centrifugal mass shed limit (rotation rate at which inwardly directed gravitational forces 
are balanced by outwardly directed centrifugal forces; £!* = fhns) (Bhattacharyya ct al. 
2000). 



EOS label 


compact star 


EOS model 


A 


SS 


Dey et al. (1998), their model SSI 


B 


SS 


Farhi and Jane (1984), B = 90 MeV/fm 3 , m s = 


C 


ss 


Farhi and Jaffe (1984), B = 60 MeV/fm 3 , m s = 200 MeV 


D 


ss 


Farhi and JafTe (1984), B = 60 MeV/fm 3 , m s = 


E 


NS 


Pandharipande (1971), hyperonic matter 


F 


NS 


Baldo et al. (1997), nuclear matter 


G 


NS 


Sahu et al. (1993), nuclear matter 



Table 1. The list of EOS models used in this work. 
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2.2. Disc temperature profile calculation 

The effective temperature of a geometrically thin blackbody disc is given by 
T cS {r) = (F(r)/a)^ (2) 
where a is the Stefan-Boltzmann constant and F, the X-ray energy flux (due to viscous 
dissipation) per unit surface area. We calculate F using the expression (Page & Thorne 
f974; valid for a geometrically thin non-self-gravitating disc embedded in a general 
axisymmetric space-time of a rotating black hole): 

m - £m (3) 

where 

f(r) = -n K , r (E - n K i)~ 2 I (E-Cl K l)l ir dr (4) 

J Tin 

Here r; n is the disc inner edge radius and a comma followed by a variable as subscript to 
a quantity, represents a derivative of the quantity with respect to the variable. 

Although formulated for the case of black holes, these expressions also hold for NS 
(Bhattacharyya et al. 2000). Depending on the mass of a compact star and the EOS 
describing it, there exist two possibilities for the location of ISCO vis-a-vis the surface 
of the star: R > r or b, or R < r or b- We assume that if R > r or b, the inner-edge of the 
accretion disc touches the star and we take r- m — R; otherwise, r; n = r or b- For strange 
stars, the inner-edge of the disc rarely touches the surface of the star (as described in 
the results section). So the above expressions will almost always be exactly valid for 
SSs. Using these results, it is straightforward to calculate the temperature profile of the 
accretion disc as a function of M and fi* of the central star for any chosen EOS. 

2. 3. Equation of state 

For strange quark matter we use two phcnomcnological models for the EOS. The first 
one is a simple EOS (Farhi and Jaffe 1984) based on the MIT bag model for hadrons. We 
begin with the case of massless, non-interacting (i.e. QCD structure constant a c = 0) 
quarks and with a bag constant B = 60 MeV/fm 3 (hereafter EOS D). Next, we consider 
a finite value for the mass of the strange quark within the same MIT bag model EOS. 
We take m s = 200 MeV and m u = m d = 0, B = 60 MeV/fm 3 , and a c = (EOS C). 
To investigate the effect of the bag constant, we take (almost) the largest possible value 
of B for which SQM is still the ground state of strongly interacting matter, according 
to the strange matter hypothesis. For massless non-interacting quarks this gives B = 90 
MeV/fm 3 (EOS B). The second model for SQM is the EOS given by Dey et al. (1998), 
which is based on a different quark model than the MIT bag model. This EOS has 
asymptotic freedom built in, shows confinement at zero baryon density, deconfinement 
at high density, and, for an appropriate choice of the EOS parameters entering the model, 
gives absolutely stable SQM according to the strange matter hypothesis. In the model 
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by Dey et al. (1998), the quark interaction is described by a screened inter-quark vector 
potential originating from gluon exchange, and by a density-dependent scalar potential 
which restores chiral symmetry at high density (in the limit ol massless quarks). The 
density-dependent scalar potential arises from the density dependence of the in-medium 
effective quark masses M q , which are taken to depend upon the baryon number density 
riB according to M„ = m„ +310MeV x sechiv— ) , where n is the normal nuclear matter 
density, q(= u, d, s) is the flavor index, and v is a parameter. The effective quark mass 
M q (n B ) goes from its constituent masses at zero density, to its current mass m q , as n B 
goes to infinity. Here we consider a parameterization of the EOS by Dey et al. (1998), 
which corresponds to the choice v — 0.333 for the parameter entering in the effective 
quark mass, and we denote this model as EOS A. 
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Fig. 1. Logarithmic plot of pressure vs. matter density for the EOS models used here. 
The density and pressure are in units of 1.0 x 10 14 g cm" 3 and (1.0 x 10 14 ) c 2 cgs 
respectively. 

For NSs, we use three representative equations of state which span a wide range 
of stiffness. A very soft EOS is the one for hyperonic matter given by Pandharipande 
(1971) which we denote as EOS E. As a representative stiff equation of state, we use the 
EOS by Sahu, Basu & Datta (1993) (hereafter EOS G). Finally, we consider the EOS 
by Baldo, Bombaci & Burgio (1997), which is a microscopic EOS for /3-stable nuclear 
matter based on the Argonnc wl4 nucleon-nucleon interaction implemented by nuclear 
three-body forces (EOS F). The latter EOS is intermediate in stiffness with respect to 
the previous two EOS models. 

A list of the designation along with the salient features of the EOS models used here 
is provided in Table 1. 
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We also display the qualitative variations in these EOS models in a log-log plot of 
Fig. 1. The differences between SS and NS EOS are plainly evident, especially at lower 
pressures. 



14 r 




2.5 



Fig. 2. Angular speed (fi*) as a function of total angular momentum (J) for strange star. 
The curves are labelled by the nomenclature of Table 1 and are for a fixed gravitational 
mass (M = 1.4 M Q ) of the strange star. 



3. The Results 

We have calculated the structure parameters and the disc temperature profiles for rapidly 
rotating, constant gravitational mass sequences of SSs in general relativity. The results 
for SS are compared with those for NS. For illustrative purposes here, we have chosen 
the value of gravitational mass to be 1.4 M© . 

Fig. 2, depicts the variation of f2* with the total angular momentum (J) for constant 
gravitational mass and for the four SS EOS. The curves extend from the static limit to 
the mass-shed limit. The striking feature here is that, although J increases monotonically 
from slow rotation to mass-shed limit, f2* shows a non-monotonic behaviour: maximum 
value of f2» (i.e. f2™ ax ) occurs at a value of J lower than that for mass-shed limit. 
Although this seems to be a generic feature for SS EOS, f2» is always a monotonic 
function of J for constant gravitational mass NS sequences and hence constitutes an 
essential difference between SS and NS (see section 4 for discussions). Our calculations 
show that at maximum f2„, the ratio of rotational kinetic energy to total gravitational 
energy: T/W approaches the value of 0.2 (see next paragraph). It has been pointed out by 
Gourgoulhon et al., 1999 that such high values of T/W make the configurations unstable 
to triaxial instability. It can also be noticed that for stiffer EOS, the star possesses a 
higher value of J at mass shed limit (ri™ ax also occurs proportionately at higher J). 
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Fig. 3. Angular speed (Q*) as a function of the ratio of rotational kinetic energy and 
gravitational binding energy (T/W) for strange star. Curve labels have the same meaning 
as in Fig. 2. 



In order to expand upon the results of Fig. 2, we plot fi* vs. T/W for various SS 
EOS in Fig. 3. It is seen that for all SS EOS, T/W becomes greater than 0.25 at mass- 
shed limit, while for NS EOS it is usually between 0.1 and 0.14 (Cook et al., 1994). 
Interestingly, for all SS EOS, flf ax occurs at about the same value of T/W (« 0.2). 

Fig. 4 displays f2* and f2i„ (i.e. the Kcplerian angular speed of a test particle at r; n ) 
against J. The four panels arc for the four SS EOS we use. We notice the interesting 
behaviour that f2* and f2; n curves cross each other at a point near fi™ ax . For rotating 
NS configurations, since the equality r; n = R is almost always (except for very soft EOS 
models: Fig. 1, Bhattacharyya et al. 2000) achieved for rotation rates well below that at 
mass-shed limit (for M = 1.4 M Q ), always fi* < fi; n (the equality is achieved only at 
mass shed limit). On the other hand, for SS, r or b is almost always greater than R (as 
explained in the next paragraph) and when the star approaches Keplerian angular speed 
at the equator, f2* becomes greater than Oi n . 

Fig. 5 is a plot of the variation of r- ln and R with O* for four SS EOS. We see that 
the behaviour of R is monotonic from slow rotation to the mass-shed limit, even though 
that of f2* is not. As mentioned earlier for all O* from static limit upto mass shed limit, 
r; n > R for 3 SS EOS. Only for the stiffest SS EOS, that we have chosen, does the disc 
touch the star (for an intermediate value of f2»). This is distinct from the case of NS 
(see Fig. 1 of Bhattacharyya et al. 2000). The reason for such a behaviour is the non- 
monotonic variation of r or b with J for SS (contrary to the case of NS and black holes); 
this is discussed further in the next section. 

In Fig. 6, we plot the variation of r- ln with fi* for three SS EOS and two NS EOS: for 
each case, our softest EOS and our stiffest EOS have been chosen. In addition, we display 
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J (10 4B gm cm 2 s-') J ( 10 49 gm cm 8 s" 1 ) 



Fig. 4. Angular speed (fi*) of the strange star (solid curve) and the Keplerian angular 
speed (Qin) of a test particle at the inner edge of the disc (dashed curve) as functions of 
total angular momentum (J) of strange star. The curves are for a fixed gravitational mass 
(M = 1.4 Mq ) of the strange star. Different panels are for different SS EOS models. 

the corresponding results for EOS model B too. It is clear that in the r- m —£l* space, there 
exists a region that is spanned by both NS and SS configurations. Interestingly, however, 
there also exists certain regions occupied exclusively by either SS or NS configurations. 
The possible observational consequences of this result is discussed in the next section. 

Fig. 7 displays the radial profiles of temperature: (i) assuming a purely Newtonian 
accretion disc and (ii) considering general relativistic accretion discs for (a) SS (EOS 
D) and (b) NS (EOS F), each represented by two configurations: the non-rotating and 
mass shed for M = 1.4 Mq . We also display the temperature profile (curve 5) for a SS 
configuration of M = 1 Mq , described by EOS (A) (the constraints obtained by Li et 
al. 1999a; 1999b) and having a period P = 2.75 ms (the mass and period corresponding 
to that inferred for the source 4U 1728-34: Mendez & Van der Klis 1999). It must be 
remembered that in this figure (and the next) , curve 5 represents the temperature profile 
for a different M value than the rest of the curves and is displayed in the same figure, only 
for illustrative purposes. From this figure we see that for M = 1.4 Mq , the Newtonian 
value of temperature is about 25% higher than the general relativistic value near the 
inner edge of the disc. This shows the importance of general relativity and rotation near 
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Q, (10 3 rad s~ l ) Q, (10 3 rad s~') 

Fig. 5. Disk inner edge radius (n n , solid curve) and strange star radius (R, dashed 
curve), as functions of angular speed f2* for various EOS models. The curves are for a 
fixed gravitational mass (M = 1.4 M ) of the strange star. 



22 




2 4 6 8 10 12 14 
n t (10 3 rad s" 1 ) 

Fig. 6. Disk inner edge radius (n n ) as a function of angular speed fi* of the compact 
star. The curves have their usual meaning. 



the surface of the star. The difference between the effects of SS EOS and NS EOS on 
temperature profiles (at the inner portion of the disc) is also prominent at mass-shed limit 
(due to the difference in rotation rates for these two configurations). Such differences in 
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Fig. 7. Accretion disc temperature profiles: Curve (1) corresponds to the Newtonian 
case, curve (2) to the Schwarzschild case (coincident curves for NS EOS model F and SS 
EOS model D), curve (3) to a neutron star (EOS model F) rotating at the centrifugal 
mass-shed limit and curve (4) to a strange star (EOS model D) rotating at the centrifugal 
mass-shed limit. For curve (1) it is assumed that, n n = 6GM/c 2 . The curves (1-4) are for 
a fixed gravitational mass (M = 1.4 M Q ) of the compact star. Curve (5) corresponds to a 
configuration that has M = 1 M Q and 17* corresponding to a period P — 2.75 ms (inferred 
for 4U 1728-34; see text) and described by EOS model A. In this and all subsequent 
figures, the temperature is expressed in units of 10 5 K, where Myj is the steady 

state mass accretion rate in units of 10 17 g s _1 . 



temperature profiles are also expected to show up in the calculations of spectra at higher 
energies. 

In the panel (a) of Fig. 8, we display the temperature profiles for configurations (as in 
Fig. 7) composed of SS EOS (D) (curves 1-4), represented by different 0* (corresponding 
to f2* = 0, minimum n n , f2* = f2™ ax and mass-shed limit); curve (5) is the same as in 
Fig. 7. The behaviour of temperature profiles is non-monotonic with il*. The panel (b) 
shows the temperature profiles at mass-shed for various SS EOS along with curve (5). 
Here the temperature profiles show monotonic behaviour with the stiffness of EOS. The 
behaviour of the temperature profiles in both the panels are similar to those calculated 
for NS (Bhattacharyya et al. 2000). Notice the substantial difference in the maximum 
temperature; sufficiently sensitive observations are, therefore, expected to complement 
the findings of Li et al. (1999a; 1999b). 

The variations of -Ed, E-ql, the ratio -Ebl/^d and T^ ax with f2* are displayed in 
Fig. 9. Each plot contains curves corresponding to all the SS EOS models considered 
here. The behaviour of all the curves are similar to those for any NS EOS (see Fig. 5 
of Bhattacharyya et al. 2000). The only difference being that due to the non-monotonic 
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Fig. 8. Temperature profiles incorporating the effects of rotation of the strange star. 
The plots correspond to (a) EOS model D and an assumed strange star mass of M = 
1.4 Mq (curves 1-4) for rotation rates: 17* = (curve 1), fi» = 3.891 x 10 3 rad s _1 (curve 
2), 0, = 7.373 x 10 3 rad s" 1 (curve 3), 0* = 7.163 x 10 3 rad s" 1 = fi ms (curve 4), (b) 
the same assumed mass and 0* = fl ms for the four EOS models (A):curve 1, (B):curve 
2, (C):curve 3 and (D):curve 4. In both panels, curve (5) is the same as that in Fig. 7. 

behaviour of from slow rotation to mass-shed limit for SS EOS, making the curves 
turn inward at the terminal (mass shed) rotation rate. 

In Fig. 10, we make a comparison between SS and NS for the same quantities displayed 
in Fig. 9. We have used three SS EOS and two NS EOS models (the softest and the 
stiffest for each case). In all the panels, SS and NS both are seen to have their own 
exclusive regions in the high and low f2* parameter space respectively. This is especially 
prominent for Sbl and Ebl/Eb- We also notice that for SS, at SI* = £l™ ax , the values 
of .Ebl ~ 0.05 and Eql/Ev s=a 1.0 for all EOS. On the contrary, for neutron stars, both 
E BL and E BL /E D become w at ft™ ax (= fi ms ). 

4. Summary and Discussion 

In this paper we have calculated the structure parameters and the disc temperature 
profiles for rapidly rotating SSs (for constant gravitational mass sequence with M — 
1.4 Mq ) and compared them with those for NSs with the aim of finding possible ways 
to distinguish between the two. For the sake of completeness, we have compared the 
properties of these two types of stars all the way from slow rotation to mass-shed limit. 

The striking feature of SSs is the non-monotonic behaviour of il* with J such that 
Sl™ ax occurs at lower value of J than that of the mass-shed limit. Hence the other SS 
structure parameters become non-monotonic functions of f2*. This behaviour is observed 
even for the constant rest mass sequences of SS (e.g. Gourgoulhon ct al 1999; Bombaci 
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Fig. 9. The variations of the En, Ebl, Ebl/Ed and T^ ax with 0* for a chosen strange 
star mass value of 1.4 M Q for the four SS EOS models. The curves have the same 
significance as Fig. 3. 



et al. 2000). In contrast, for NSs, the structure parameters are all monotonic functions 
of f2*. An implication of the non-monotonic behaviour of fi* with J is that if an isolated 
sub-millisecond pulsar is observed to be spinning up, it is likely to be a SS rather than 
aNS. 

Because of higher values of T/W (£0.2), SSs are more prone to secular instabilities 
compared to NSs at rapid rotation (Gourgoulhon et al. 1999). Our calculations show that 
at f^ lax , T/W > 0.2. 

Another important feature of SS gravitational mass sequence (in contrast to the 
corresponding NS sequences) is the crossing point in ft* and ft; n . This feature has im- 
portant implication in models of kHz QPOs: for example if ft* is greater than Sl- m , the 
beat-frequency models ascribing higher frequency to Keplerian frequencies will not be 
viable. 

It can be noted from Fig. 4, that with the increase in stiffness of the EOS models, J cro ss 
increases and ft* lCr0 ss (the subscript "cross" corresponds to the point ft; n = ft*) decreases 
monotonically. It is also seen that in general all the quantities vary monotonically with 
the stiffness for both SS and NS EOS (see also Bhattacharyya et al. 2000). 
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Fig. 10. Same as Fig. 9, except the fact that here two NS EOS models and two SS EOS 
models are used. The curves have the same significance as Fig. 6. 



For SSs, the inner-edge of the accretion disc rarely touches the surface of the star 
(even for maximum rotation rates), while for rapidly rotating NSs, the accretion disc 
extends upto the stellar surface for almost all rotation rates. Since the inner accretion 
disc boundary condition is different for both these cases, we expect important observable 
differences (both temporal and spectral) in X-ray emission (from the boundary layer and 
the inner accretion disc) from SSs and NSs. 

A brief note on the variation of r orb , with specific angular momentum is in order here. 
As mentioned earlier, beyond a certain value of the angular momentum, the radius of 
the ISCO increases with increasing angular momentum - a property not seen either in 
the case of NSs or black holes. The reason for this can be traced to the radial gradient 
of the angular velocity of the particles at the marginally stable orbit and the analysis is 
described as follows: 

For the metric described by equation (1), the second derivative of the effective po- 
tential may be written as (Cook et al 1994; Thampan & Datta 1998) 

V X r = 2 



\^r-l 2 r)-\^ 2 y-\r + - r 



v 2 



+ [2 + r(7 >r — X tr )]vv. r — e uj r rv 

+^(7 2 r -^)-^V r = (5) 
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Simplification of this, using the other equations of motion (Bardeen 1970), yield 



r 2 (l - v 2 )V rr = -X 



2 



rfl r 1 — V' , 



(6) 



Where X = v 2 (2 + r7. r — rp /r ) + r(j <r + p, r )- The marginal stability criterion, therefore, 
yields the rate of change of the marginal stable orbit, with respect to specific angular 
momentum (J = J/M 2 ) as: 

r orb j = r orh +2^— A\ (7) 

[ n - uj ri jr v(v 2 - 1) x j 

where the terms in the parenthesis are to be evaluated at r or b- We calculate the four 
terms in the parenthesis in equation (Q) and find that the second term dominates the net 
rate of change of r or b with j. Which implies that at the value of j where r or bj changes 
sign, although the first three terms are observed to change sign, the net sign is only 
dependent on that of il. r j at ISCO. 

From Fig. 6 we see that for f2* in the range (0,4028) rad. s" 1 (the second quantity in 
the range is the rotation rate of PSR 1937+21: Backer et al. 1982, the fastest rotating 
pulsar observed so far), a major portion of the ri n -0» space is occupied exclusively by 
NS. So if r; n can be determined independently from observations (for example, by fitting 
the soft component of the observed spectrum by the XSPEC model "disebb" available 
in XANADU: see for example Kubota et al 1998, or, from the observed kHz QPO fre- 
quencies), there is a fair chance of inferring the central accretor to be a NS rather than 
SS (provided the mass of the central accretor is known by other means). This is also 
applicable to Ebl and Ebl/Ed (Fig. 10). It is also to be noted that Li et al. (1999 a; 
1999b) did a similar search in the M — R parameter space and concluded the millisecond 
X-ray pulsar SAX J1808.4-3658 and the central accretor in 4U 1728-34 to be likely SSs. 
If, indeed this is true, then it is possible to constrain the stiffness of the equation of state 
of SQM (Bombaci 2000), and to exclude EOS models (like EOS C and EOS D) stiffcr 
than our EOS B. 

Calculation of the accretion disc spectrum involves the temperature profiles as inputs. 
The spectra of accretion discs, incorporating the full effects of general relativity for NSs 
(Bhattacharyya et al. 2000;2001a;2001b) show sensitive dependence on the EOS of high 
density matter. The similarity in the values of the maximum disc temperature implies a 
similar indistinguishability in the spectra of SS as compared to NS. Nevertheless, just as 
Ebl and other quantities show that NS exclusively occupy certain regions in the relevant 
parameter space, we expect that it will be possible to make a differentiation between 
these two compact objects by modelling the boundary layer emission. If as mentioned 
in previous paragraph, we exclude EOS models stiffer than B, then from Fig. 10, we see 
that a fairly accurate measurement of Ebl (Lbl/ Mc 2 ) and f2* can indicate whether the 
central accretor is a NS or SS if the corresponding point falls outside the strip defined 
by curves B and E. 
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The current uncertainties in theoretical models of boundary layer emission and the 
variety of cases presented by models of rotating compact objects, calls into order, a 
detailed investigation into these aspects of LMXBs - especially with the launch of new 
generation X-ray satellites (having better sensitivities and larger collecting areas) on the 
anvil. 
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